function answer = fact(n)
      
    if (n == 0)
	 answer = 1;
    else 
	 answer  =  n * fact(n-1);
    endif

endfunction